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Abstract. - We show that the low-temperature physics of molecular nanomagnets, contrary 
to the prevailing one-molecule picture, must be determined by the long-range magnetic order- 
ing due to many-body dipolar interactions. The calculations here performed, using Ewald's 
summation, suggest a ferromagnetic ground state with a Curie temperature of about 130 mK. 
The energy of this state is quite close to those of an antiferromagnetic state and to a glass of 
frozen spin chains. The latter may be realized at finite temperature due to its high entropy. 



In recent years, molecular nanomagnets Mni2-Ac and Feg of spin 10 have been 

the object of intense theoretical and experimental studies due to their interesting magnetic 
properties. They exhibit spectacular quantum effects, such as regular steps in the mag- 
netic hysteresis [^|^], quantum interference effects in the magnetic relaxation Q, quantum 
avalanches |^,|^ , coherent quantum oscillations of spin between two classical energy minima |^ , 
first- and second-order crossover between thermal and quantum regimes Q, etc. Most of these 
effects can be observed in macroscopic measurements and can be studied by quasi-classical 
methods, which places the physics of high-spin molecular clusters in the domain of macro- 
scopic quantum tunneling [^^[ll]] . A significant theoretical effort has been mounted to explain 
the observed phenomena [[l2| . Most of the theoretical works to date treated Mni2-Ac and Fcg 
crystals as consisting of non-interacting spin-10 clusters, addressing basically the question 
of transitions between spin levels of individual clusters due to internal and external fields, 



phonons, and nuclear spins. More recently Prokof'ev and Stamp |13| suggested that in cer- 
tain experiments the dipolar interactions between clusters can be responsible for the relaxation 
law (see also |lj, ^). This suggestion has been confirmed by computer simulations [p^p^ 
and by the experimental evidence of memory effects 0,11^ recently observed in both Mni2-Ac 
and Fes. 
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The question addressed by us is this: can dipolar interactions induce magnetic ordering in 
Mni2-Ac and Fes? If this was true, the low temperature physics of molecular nanomagnets 
would be the physics of magnetic order in a many-body system, which would be more inter- 
esting than the current spin-10 physics. It is important to emphasize that the possibility of 
ordering has been overlooked in previous works studying, e.g., the relaxation of these systems. 
For instance, the hole-digging method can be regarded, in the light of this paper, as due to 
irreversibility and memory effects in a spin-glass phase (a well-known effect |19|). The subject 
of magnetic ordering induced solely by dipolar interactions is not newpQ-p^. Experiments 
established magnetic order in dilute rare earth insulating compounds [p4[. Theoretical studies 
performed to date were restricted mainly to cubic lattices. The main theoretical result is 
that the lowest energy configuration for a simple cubic (sc) lattice corresponds to an antifer- 
romagnetic configuration (with ordering vector k = (0.5,0.5,0)), while in the body-centered 
(bcc) and the face-centered (fee) cubic case, the ferromagnetic configuration is energetically 
favorable. 

Both Mni2-Ac and Fes have strong easy-axis anisotropy. The H = problem is then the 
one of an Ising system interacting via magnetic dipole forces. There is a significant difference 
between the two systems, however. Firstly, in Mni2-Ac the terms breaking the uniaxial 
magnetic symmetry aX H — are weak. Consequently, at low temperature (T < IK), the 
spins of Mni2-Ac molecules are blocked along the anisotropy directions chosen by the history 
of the cooling process. In addition, the hyperfine fields in Mni2-Ac are as strong as the dipole 
fields but random from one Mni2-Ac cluster to another. Consequently, a highly metastable 
spin-glass order, if any, should be expected in Mni2-Ac as temperature is lowered. 

In Feg the terms violating the uniaxial symmetry are large, leading to thousands of quan- 
tum transitions per second between the equivalent easy-axis directions at T = |^]. Also, 
Fes has smaller lattice parameters (a = 10.522(7) A, 6 = 14.05(1) A, c = 15.00(1) A, 
a = 89.90(6)°, /3 = 109.65(5)°, 7 = 109.27(6)° |||), making dipole interactions more im- 
portant, while hyperfine interactions in Feg are very weak. The Fes cluster is chemically 
constituted by an octameric Fe^^^ cation, in which six of the Fe^^^ atoms coordinate with a 
cyclic amina. The cyclic amina ligands form the "surface" of the complex cation, effectively 
isolating the Fe^^^ atoms from the surrounding (more details on the chemical structure of Feg 
can be found in Ref. j25|) . Therefore, other interactions between Fcg clusters, such as ex- 
change or super exchange, are also negligible. Thus, in Fcg the low temperature magnetically 
ordered ground state due to dipole interactions must be simpler and much easier to achieve 
than in Mni2-Ac. The estimate of the temperature at which dipolar interactions might begin 
to be of importance can be obtained from the order of magnitude of the dipolar energy. One 
gets Tdip ~ Ed = igi^BS^^/Vc « 130 mK, where /iB is the Bohr magneton, g « 2 is the 
gyromagnetic ratio |g,^,|27| and 14 ~ 1956 A'^ is the volume of the unit cell |25| of Feg. Note 
that if the ordering is ferromagnetic (as will be shown below) , deviations of the actual critical 
temperature from this mean-field estimate should not be too large, as the upper critical di- 
mension for an Ising dipolar ferromagnet is d = 3 |p3| . In this work we shall find the classical 
ground state of Feg due to dipole interactions, assuming that the only role of the terms in the 
Hamiltonian which are responsible for tunneling between individual easy directions, is to help 
to achieve that ground state. 

To find, through direct minimization, the lowest energy configuration, one can study the 
Fourier transform of the dipolar interaction. Here we are assuming for simplicity that each 
Feg cluster is characterized by a single dipole moment located at the center of the cluster. 
Given that magnetic moments of eight iron ions belonging to an Feg cluster are separated by 
an appreciable distance, this, of course, is a crude approximation. We believe, however, that 
it catches the physics of the problem. 
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The Hamiltonian we consider for the crystal, assuming that the external field is zero, is 

'^dipolc ^" "^single — ion i 

(1) 

where ^ 

Wdipolo = ~ 2 51 ^al3,ijS°'S^ (2) 

and 

Wsi„glc-ion=^-A-(5f +65f') . (3) 



For Fc8, one has K w 230 niK, b « 0.4 and Eb « 130 mK |2||26|,g7[. We will consider the limit 
in which the spins align themselves along the easy axis (the z axis in our notation), which we 
take to be approximately parallel to the crystallographic axis a pOt - Thus we can project each 
spin onto the basis 15*2 = ±5). Note that the energy difference for a small angle rotation from 
15^ = S) to \Sz ^ S - 1) is AE w K{2S - 1) w 5 K > £'d in Fcg, so this approximation is 
quite accurate in the relevant range of temperature. Taking into account these considerations, 
we obtain an Ising-like array of spins with the effective Hamiltonian 

w=-^E^--^^--^;+E^<' (4) 

where Dzz,ij = [3cos^0y — % is the angle between and the easy x-axis and 

A is the tunnel splitting responsible of tunneling between the \Sz ~ ±5") states. Here 
is a Pauli matrix for the site i. If the main source of tunneling is the Sf'^ term, A is of 
order Kb^ [^(for Fcg a typical value is A « 10^^ K <C Ed)- This Hamiltonian describes the 
" quantum spin glass" problem |2^,^^ . From now on, we will restrict our considerations to the 
classical approximation, i.e., S oo. We will accordingly neglect the tunneling term |(^)| and 
regard rf as a number. Thus, the problem reduces to finding the configuration {rf } which 
minimizes 

^ = — ^ E ' (5) 

with the restriction 

rf = ±1, Vz . (6) 
We consider the Hamiltonian in the momentum space 

n/E^ = J d^fc T^{k)D.Ak)T^-k) , (7) 

where the integration is performed over the first Brillouin zone and T^{k) and Dzz{k) are the 
Fourier transforms of and Dzz,ij, respectively, given by 

i 

Dzz{k)=Y,Dzz.oe'''''^ . (9) 



(^) The tunneling term is necessary to obtain the true quantum ground state. Even at finite temperatures, it 
induces fluctuations that might destroy the order found here. 



4 



EUROPHYSICS LETTERS 



Here we have made use of the translational invariance of Dzz.ij- The sum in Eq. converges 
very slowly. In fact, it is non-analytical as fe — s- 0, this being a manifestation of the shape 
dependence of the energy in the presence of a net magnetization. This sum has been computed 
using Ewald's method [^,31,3^. In this method, the Fourier transform of the dipolar tensor 
is given by 



exp(-7r262/i?2^ 



Gi 



(3 + 2a;J - gaf) 



' 3V^ 



Ja0 



9aP 



erfc(a;i) > exp(z27rfc^rf ) 



(10) 



where Gi runs over the reciprocal lattice (the prime indicates that we exclude the term for 
Gi = (0,0,0)), ri runs over the Bravais lattice, Xi = Rri, bi = Gi + k and gap is the metric 
associated with the unit cell basis. Note that the basis of the reciprocal lattice used to expand 
the reciprocal space vectors in components is the dual basis e" of the real space unit basis ep, 
defined by e"{ep) = Sp (and not by the more conventional relation e°'{ep) = 27r5^). i? is a 
real parameter which controls the convergence of the series (the value of D is independent of 
it). For a cubic lattice, the optimal value is i? = 2.0/a, whereas for Feg it has been found that 

1/3 

convergence is fastest for both the reciprocal and the real space sum using R ~ l.8/Vc ■ 

Now it must be noted that one cannot realize spin structures for arbitrary k with Ising 
spins. Indeed, in order to satisfy restriction (^, k must be either in the center of the Brillouin 
zone (ferromagnetic arrangement) or at its boundary, i.e., we must have = 0,1/2, with 
IJ, — 1, 2, 3. Moreover, the dipolar interaction energetically favors the ferromagnetic alignment 
of spins in chains along the a axis. Thus, we will restrict our consideration to the case A:i = 
(configurations with /ci = 1/2 have been found to have much higher energy, as expected). The 
distinct configurations we need to consider are 



fco = (0,0,0), 
fc2 = (0,1/2,0), 



fci = (0,0,l/2), 
fe3- (0,1/2, 1/2) 



(11) 



Here fco corresponds to a ferromagnetic arrangement, and fci, k2 and k^ correspond to anti- 
ferromagnetic arrangements with planes of alternating magnetization ([001], [010] and [011]- 
planes, respectively). The energy per spin of these configurations can be computed using 
Eq. ( p^ ) and the relation E = —{l/2)Dzz{k)E'o. One obtains 



Eq = -4.10 £^0 
E2 ^ -Am Ed 



El : 

^3 



-4.05 £^0 
-3.97 



(12) 



Therefore, we find that the lowest energy configuration corresponds to ferromagnetic alignment 
of the spins, with an energy per spin Eq sa —530 mK (note that due to the smallness of the 
energy differences, unaccounted effects may, in principle, shift the ground state; this will be 
discussed elsewhere). In evaluating Eq, one must note that the non-analyticity of the first 
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term in Eq. ( |To| ) can be attributed to shape-dependent demagnetizing effects, as discussed 
above. We can define the demagnetizing factor as 



L = 4tt hm . (13) 

fc^(o,o,o) \ k J 

Then, one must talce L — for needle-shaped samples (or for systems with periodic boundary 
conditions without any free surfaces), which is the situation we have assumed in computing 
Eq. Nevertheless, the bulk free energy of the system should be shape independent in the 
absence of an external field (see, e.g., Refs. |^,|3^]). Physically, for other sample shapes the 
demagnetizing energy is reduced due to the formation of domains in the ferromagnetic state. 
In this case, though, we should include a domain wall energy contribution. However, this 
contribution can be shown to be negligible in the thermodynamic limit due to the macroscopic 
size of the domains [|5|, and hence Eq should be shape independent. 

Besides these configurations, one must note that no general arguments exist which would 
rule out the existence of lower lying ground-states with structures having larger unit cells (such 
as arrangements where two planes of ferromagnetically aligned up-spins alternate with two 
planes of down spins). However, numerical investigation of Dzz{k) by standard optimization 
methods (e.g., amoeba method), seems to indicate that k = ko is the global maximum of 
Dzz- As any configuration can be expanded in the "plane wave" fc-basis, we conclude that 
the ferromagnetic arrangement is, effectively, the lowest energy configuration. 

Particularly interesting and physically relevant (because of its large entropy) is also a "spin 
glass" -like phase, where chains of ferromagnetically aligned up-spins in the a-direction alter- 
nate with chains of down-spins, such that the spin arrangement in the 6c-plane is completely 
random. The energy of this phase is given by 

Eso/E^^-l ^ = -2C(3)^--4.04, (14) 

where C(p) = y^^-, n~P is the Riemann zeta-function. 

In recent work |Q , it has been claimed that the easy axis is parallel to the crystallographic 
b axis. For the sake of completeness and in order to compare the results here obtained with 
those presented in Ref. [^6| , we have considered alternative orientations of the easy axis. 
Taking the easy axis parallel to the crystallographic b axis, one obtains that the lowest energy 



configuration corresponds to ki — (0,0,1/2), in agreement with MC simulation |36|. The 
energy here computed is Ei/E-o ~ —1.78, which differs by less than 10% from the value quoted 
in Ref. js^ . This small difference could be attributed to finite-size effects and the truncation 
of the interaction used in MC simulations. In the case that the easy axis is taken parallel to 
the c axis, the ground state configuration is k^ = (1/2, 1/2,0), with energy E^/Ed = —1.55. 

Summarizing, the present work suggests the onset of the ferromagnetic ordering in Fes 
due to dipolar interactions below T « 130 mK. The effect should exist in magnetic fields 
below the dipole field for which we obtain H = D zz{ko)g hbS /Vc ~ 77 mT. It can be observed 
by conventional methods. The most direct method would be the observation of a different 
relaxation dynamics above and below Tc for the sample initially magnetized along the mag- 
netic anisotropy axis. After the field is turned off, such a sample should relax towards lower 
magnetization above Tc and towards higher magnetization below Tc. Measurements of the 
specific heat of Feg down to very low temperatures would be helpful as well. However, achiev- 
ing the ferromagnetic phase might be difficult, due to its low entropy and the presence of high 
anisotropy barriers, and experiments might probe a spin-glass ordering, instead. In that case 
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one should look for memory effects. In fact, aging and rejuvenation phenomena observed at 
mK temperatures in Mni2-Ac and Fes (l^'llll very similar to those found in spin-glasses 
and disordered ferromagnets [ p^ . Another interesting problem, not addressed by this work, is 
the spectrum of collective excitations in the magnetically ordered phase of Fes ■ The presence 
of the tunneling term in the Hamiltonian, suggests the possibility of a spin-flip excitation 
( "tunnelon" ) . This and other questions will be studied elsewhere. 

We are grateful to Wolfgang Wernsdorfer for discussion and for bringing our attention to 
Ref. . XMH acknowledges support from the Generalitat de Catalunya. EMC acknowledges 
NSF grant No. 9978882. This work had begun during the authors' stay at the Institute of 
Theoretical Physics in Santa Barbara, within the Program on Magnetic Phenomena in Novel 
Materials and Geometries. 
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